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FOREWORD 


This  report  was  prepared  for  the  Air  Force  Flight  Dynamics 
Laboratory,  Flight  Research  and  Test  Branch,  Air  Force  Systems  Command, 
by  the  Flight  Systems  Division,  The  Bendix  Corporation,  Teterboro, 

New  Jersey  under  Air  Force  Contract  No.  F33615-72-C-1753,  data  document 
item  A006,  Improved  Flare  Law  Study.  This  work  was  performed  in  conjunction 
with  other  program  tasks  during  the  time  period  January,  1972,  through 
December,  1974.  The  Air  Force  Program  Manager  was  Capt  T.  Imrich  with 
Project  Engineers,  Lt  R.  P.  Denaro  and  Lt  B.  Kunclw. 

The  Flight  Systems  Division  effort  was  under  the  direction  of 
Mr.  F.  G.  Adams,  principal  investigator.  The  simulation  and  analysis  task 
was  conducted  by  Mr.  V.  Muehter  with  contributions  by  Mr.  K.  Moses, 

Assistant  Chief  Engineer.  Supporting  data  and  configuration  requirements 
were  supplied  by  Messrs.  J.  Woloshen,  M.  Sforza,  and  S.  Skaritka. 
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SECTION  1 
INTRODUCTION 


This  study,  in  response  to  data  document  item  A006  was  per- 
formed to  investigate  improved  control  law  concepts  in  the  flare  portion 
of  the  automatic  approach. 

The  present  exponential  flare  law,  using  altitude  ( h)  and  de- 
rived altitude  rate  ( H)  yields  satisfactory  results  under  most  conditions. 
It  was  felt  that  in  using  information  now  available,  such  as  range  (R)  and 
flight  path  angle  ( y),  reductions  could  be  made  in  touchdown  dispersions 
presently  occurring  under  certain  wind  shear  conditions. 

A study  was  also  performed  in  the  area  of  advanced  throttle 
control  as  a modification  to  the  present  open  loop  throttle  retard  during 
flare.  This  effort  was  directed  toward  optimizing  airspeed  bleeds  during 
flares  under  varying  wind  shear  conditions. 
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SECTION  2 


DISCUSSION 


The  present  KC-135  automatic  flare  configuration  uses,  h + 
tH  = 0 as  its  control  law.  This  forces  an  exponential  altitude  characteris- 
tic with  respect  to  time.  It  is  noted,  however,  that  no  restoring  term 
exists  should  the  aircraft  deviate  from  the  nominal  flare  path  in  space,  due 
to  gusting  or  windshears.  The  control  law  merely  requires  a particular 
sink  rate  at  each  altitude.  Thus  consistent  sink  rates  at  touchdown  are  its 
main  characteristic  with  a considerable  sacrifice  in  touchdown  dispersions 
under  varying  wind  conditions. 

1 . R,  y Flare  Law 

The  flare  law  herein  investigated,  as  an  alternative  to  the 
present  system,  is  depicted  in  Figure  1.  This  is  basically  a y control 
system  which  indirectly  controls  sink  rate  at  touchdown  with  small  vari- 
ations from  nominal,  corresponding  to  changes  in  ground  speed  due  to 
wind  conditions.  These  variations,  in  the  order  of  14%  for  a 20  knot  wind, 
are  considered  a small  price  to  pay  for  a system  which  does  not  deviate 
from  the  nominal  flare  path  to  maintain  a particular  sink  rate  during 
wind  shears.  The  y command  term  can  be  considered  to  be  a direct 
function  of  altitude  when  the  aircraft  is  tracking  the  nominal  flare  path 
( hj  = 0)  , The  h command  term  ( fi)  is  a computed  function  of  range 
(referenced  to  a point  on  runway  1000  feet  behind  the  glide  slope  trans- 
mitter) corresponding  to  the  nominal  no  wind  flare  path  in  space.  This 
additional  computed  h error  ( ht  ) term  prevents  the  system  from  allow- 
ing the  aircraft  to  fly  a path  parallel  to  the  nominal  flare  path,  should  it 
be  displaced  therefrom  for  any  reason  ( see  Appendix) . 

To  investigate  the  feasibility  of  such  a system,  comparisons 
were  made  between  it  and  the  present  exponential  system  under  various 
wind  conditions.  The  wind  models  used  were  those  suggested  by  the  Air 
Force  Flight  Dynamics  Laboratory  and  consisted  of  constant  head  and  tail 
winds,  linear  shears,  logarithmic  shears  and  knife-edge  shears.  The 
flare  responses  for  the  nominal  and  proposed  configuration  are  shown  in 
Figures  2 and  3,  respectively,  for  the  20  knot  headwind,  no  wind,  and  20 
knot  tailwind  cases.  As  expected,  the  sink  rate  at  touchdown  for  the  ex- 
ponential flare  was  held  constant  at  the  expense  of  longitudinal  disper- 
sions. The  proposed  R,  y flare,  on  the  other  hand,  allowed  approximately 
a 15%  variation  from  its  nominal  (no  wind)  sink  rate  but  displayed  very 
little  dispersion.  Figure  4 displays  the  flare  performance  in  20  knot  head 
and  tail  winds  sheared  out  linearly  from  500-0  feet.  It  is  noted  that  the  R, 
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■y  flare  has  displayed  considerably  less  touchdown  dispersion  and  at  the 
same  time  resulted  in  less  variation  in  touchdown  sink  rate.  It  was  found 
that  in  cases  where  winds  were  sheared  to  neai  zero  at  ground  level,  the 
R,  V flare  configuration  not  only  yielded  smaller  longitudinal  dispersions 
but  also  less  variations  in  sink  rate  than  the  exponential  system.  The  lat- 
ter is  due  to  the  combination  of  the  aircraft'  s increased  ability  to  remain 
on  the  nominal  y path  and  the  absence  of  wind  at  ground  level.  These  same 
characteristics  are  demonstrated  for  the  case  of  20  knot  winds  logarith- 
mically sheared  out  from  500-0  feet  in  J'igure  5. 

Knife-edge  shears  at  1 1 0 feet  from  20  knots  shearing  to  13  knots 
were  simulated  during  flares  in  Figure  6.  Again  touchdown  dispersion 
characteristics  for  the  R,  y flare  are  far  superior  to  those  of  the  exponen- 
tial flare.  This  type  of  shear  was  found  to  cause  the  greatest  variation  in 
touchdown  sink  rate  for  the  R,  y system.  Even  in  this  case,  however,  the 
choice  to  accept  a somewhat  larger  sink  rate  variation  to  ensure  far  smal- 
ler longitudinal  dispersions  is  clear. 

The  results  of  all  the  above  referenced  simulated  flares  are  sum- 
marized in  the  table  of  Figure  7.  The  overall  superiority  of  the  proposed 
R,  y flare  system  is  quite  evident. 

2.  Throttle  Retard  Considerations 

In  the  present  KC-135  throttle  system,  the  autothrottle  output  is 
frozen  at  zero  at  the  initiation  of  flare.  At  this  time  an  open  loop  ramp  is 
fed  to  the  throttle  servo  as  a constant  rate  throttle  retard,  up  to  a fixed 
limit.  The  rate  and  limit  have  been  chosen  to  provide  an  optimum  airspeed 
bleed  during  flare,  under  the  nominal,  no  wind  condition.  An  optimum 
bleed  in  this  case  i s'approximately  15  knots,  which  would  still  allow  a suf- 
ficient margin  above  the  stall  airspeed.  Referring  back  to  Figure  2,  it  has 
been  shown  that  for  20  knot  head  and  tail  winds,  the  airspeed  bleed  can  vary 
from  1 9 to  12  knots  respectively.  This  is  due  to  the  fact  that  the  time  for 
flare  is  larger  in  a headwind  than  tailwind  and  the  rate  of  bleed  is  the  same 
in  each  case. 

These  airspeed  bleed  variations  would  suggest  the  use  of  some 
form  of  closed  loop  around  the  throttles  based  upon  a comparison  between 
^round  speed  ( Vg)  and  airspeed  ( Va)  . Studies  have  shown,  however,  that 
the  slow  response  of  airspeed  to  throttle  position  precludes  the  use  of  a 
closed  loop  linear  system.  Therefore,  various  non-linear  schemes,  in 
which  the  retard  rate  and  limit  are  functions  of  ( Va-Vg),  were  investigated. 

One  configuration  intended  to  minimize  the  airspeed  variations  of 
Figure  2,  is  shown  in  Figure  8.  The  presence  of  either  a head  or  tail  wind 
is  sensed  by  comparing  Vg  with  Va  and  the  retard  limit  is  adjusted  accord- 
ingly. The  results  of  the  present  system  flares  incorporating  this  feature 
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DGE  SHEAR 


FIGURE  7.  FLARE  PERFORMANCE  SUMMARY 
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I 10' 

-1650 

-2.  0 

-830 
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♦Distance  Referenced  to  Glide  Slope  Transmitter 
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are  displayed  in  Figure  9 under  various  wind  conditions.  It  is  noted,  the 
characteristic  of  more  airspeed  bleed  in  a headwind  has  still  remained 
despite  the  fact  that  the  throttles  have  not  retarded  as  far.  This  is  caused 
by  the  fact  that  the  time  for  flare  has  been  affected  by  the  level  of  throttle 
retard,  through  the  pitching  moment  due  to  thrust  ( . Figure  10  de- 

scribes another  method  controlling  throttle  retard  by  adjusting  retard  rate 
as  a function  of  ( Va-Vg) . Its  responses  are  shown  in  Figure  11  and  show 
similar  characteristics  to  those  of  Figure  9.  Therefore,  it  seems  that  any 
attempt  to  alter  airspeed  bleed  to  compensate  for  a constant  wind  condition, 
by  altering  throttle  retard  is  self  defeating. 

Compensation  for  wind  shears,  during  flares,  may  still  be  possi- 
ble. This  has  been  an  area  of  prime  concern.  The  occurrence  of  a knife- 
edge  shearing  out  of  a headwind  or  shearing  in  of  a tailwind  may  reduce  air- 
speed close  to  stall  levels  using  an  open  loop  throttle  retard.  Therefore,  the 
configuration  of  Figure  12  has  been  under  study.  The  object  of  this  scheme 
is  to  r^-duce  the  nominal  retard  limit  as  a function  of  only  sudden  drops  in 
airspeed.  These  would  be  caused  by  airspeed  changes  too  rapid  to  allow  the 
aircraft  to  maintain  its  airspeed,  thus  the  use  of  the  washout.  Therefore, 
the  time  of  flare  would  not  be  affected,  as  discussed  above,  unless  a serious 
wind  shear  such  as  a knife-edge  were  to  occur.  In  that  event  the  limit  would 
be  reduced  to  help  maintain  airspeed.  The  responses  of  this  configuration 
and  that  of  the  present  system,  to  20  knots  headwind  with  10  knots  knife-edge 
sheared  out  at  30  feet,  are  shown  in  Figure  13.  It  is  noted  that  improvements 
can  be  seen  due  to  the  reduced  retard  limit.  As  the  shear  is  encountered,  the 
retarded  throttles  are  driven  back  toward  positive  to  oppose  the  instantaneous 
airspeed  loss.  This  results  in  a total  airspeed  loss  comparable  to  that  of  the 
nominal  no  wind  flare.  Sink  rate  and  touchdown  dispersion  have  also  been 
significantly  improved.  This  shear  occurred  at  a point  approximately  mid- 
point during  the  flare.  It  is  noted  that  shears  occurring  closer  to  the  ground 
could  not  be  compensated  for  as  well  since  the  resulting  forward  throttle  mo- 
tion would  have  less  time  to  effect  the  overall  problem. 
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VARIABLE  RETARD  LIMIT  RESPONSES-CONFIGURATI 
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FIGURE  11.  VARIABLE  RETARD  RATE  RESPONSES 
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FIGURE  13.  FLARE  RESPONSES-KNIFE -E DGE  SHEAR 


SECTION  3 


SIMULATION 


The  followiijg  perturbation  equations  were  used  to  simulate  the 
three  degree  of  freedom  longitudinal  axis. 


0 = M^w  + M^w  + M^e  + M^u  + 6T  + 6e  + H + 6s 


w = z w + z.w  + (V+z)0+zu  + z,  6e  + z H 
w w ' q'  u 6e  H 

i = g0tx..,wtx^u*Xj^6TtXj,^6e  + XjjH 
Parameters  are  defined  in  Figure  14. 

The  ground  effect  data  obtained  from  Wright  Field  was  as 

follows: 


-0.138  (1-A/A^) 

+ 0.  127  ( 1 - A/A^) 

G 

-0.  045  ( 1 - A/A) 

G 

wnere  is  the  ettective  aspect  ratio  in  the  presence  of  the  ground.  The 
quantity,  1-A/A  , designated  by  H,  is  shown  in  Figure  15,  as  a function 
of  C.  G.  altitude . 

The  horizontal  stabilizer,  denoted  above  by  6 , was  simulated 
with  a time  rate  of  0.  15  degrees  of  surface  per  second.  ® It  was  engaged 
whenever  elevator  deflection  exceeded  1.5  degrees.  The  effectiveness 
of  the  stabilizer  in  producing  a pitching  moment  was  taken  to  be  three 
times  that  of  the  elevator. 

The  twelve  flight  conditions  at  which  existing  aero  data  was  de- 
rived are  not  representative  of  the  gross  weight  to  airspeed  ratio  that  is 
currently  being  flown  on  approach.  For  this  reason,  a new  flight  condi- 
tion was  developed  with  a gross  weight  of  160,  000  pounds  and  airspeed  of 
155  knots.  Figure  14  shows  the  dimensional  stability  derivatives  obtained 
from  Wright  Field  for  Flight  Condition  2 defined  with  the  same  gross 
weight  mentioned  above,  but  an  airspeed  of  140  knots.  The  values  for 
derivatives  in  the  new  flight  condition  were  obtained  by  multiplying  by  the 
appropriate  factor,  155/140  or  (155/140)2  as  shown  in  Figure  14.  For 
uniformity,  this  flight  condition  was  simulated  throughout  the  longitudinal 
study. 
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FIGURE  14.  LONGITUDINAL  DIMENSIONAL  STABILITY  DERIVATIVES 


Parameter 

Definition 

Value  at  Flight 
Condition  2 

Value  this 
Flight  Condition 

M 

w 

ps  Vc 

m_ 

a 

-. 6219  X lO'^ 

-.  6916  X 10  ^ 

M\v 

2 

psc 

4Iy 

-.  1466  X 10  ^ 

-.  1466  X 10  ^ 

M 

q 

ps  Vc 

- ^rnq 

8554 

-.  9512 

Mu 

ps  Vc 

ly 

0 

0 

MgT 

. 1852  X 10*^ 

. 1852  X lo'^ 

>^6e 

2 

ps  V c 

. ^m5e 

1.  124 

1.  390 

Mh 

ps  V^c 

2Iy 

-. 2576 

Zw 

-ps  V 

( Clq.+Cj^) 

-.  637 

-.  708 

2m 

Zw 

-psC 

4m 

^Ld 

-.0105 

-.  0105 

-ps  Vc 
4m 

CLq 

-6.  15 

-6.  84 

I9 


Parameter 

Definition 

Value  at  Flight 
Condition  2 

Value  this 
Flight  Condition 

Zu 

m ^ 

-.  272 

302 

^6e 

-psV^  C 
2m  6e 

7,46 

9.  23 

-psV^  AC, 
2m 

-5.  113 

. 064  7 

. 0719 

Xu 

lfis_V  (VC  ) 
m 

0426 

0474 

^6T 

J_ 

m 

. 2013  X lo"^ 

. 2013  X 10*^ 

^6e 

^ C 
2m  Se 

1,  49 

1.  84 

Xh 

AC 

2m 

1.812 

08 

3Ms^ 

3.  372 

4.  170 
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The  characteristics  of  the  logarithmic  shears  used  in  the  above 
analyses  are  shown  in  Figure  16. 


SECTION  4 


CONCLUSION 


This  study  has  investigated  the  feasibility  of  improving  the 
flare  performance,  during  automatic  approach,  through  the  use  of  in- 
ertial flight  path  angle  (y)  and  range.  A comparison  was  made  between 
it  and  the  present  system,  which  uses  an  exponential,  altitude,  altitude 
rate  control  law. 

The  proposed  system  which  commands  a y,  based  upon  altitude 
and  range  information,  was  found  to  yield  improved  performance  under 
varying  wind  conditions.  Longitudinal  dispersions  at  touchdown  were 
found  to  be  much  reduced  in  all  cases  with  a slight  improvement  in  over- 
all touchdown  sink  rate  variation. 

A study  was  also  performed  in  an  effort  to  optimize  the  air- 
speed bleed  during  flare,  for  varying  wind  conditions,  through  the  use  of 
range  rate  information.  Having  investigated  various  forms  of  closed  loop 
throttle  retard  systems,  it  was  concluded  that  a non-linear  system  which 
varies  the  throttle  retard  limit  as  a function  of  ground  speed  and  airspeed 
is  the  most  promising.  However,  it  must  not  affect  throttles  on  a long  term 
basis  since  flare  performance  would  then  become  a function  of  winds  due  to 
pitching  moment  due  to  thrust  ( . Therefore,  a system  adjusting  re- 

tard limit  on  a short  term  basis,  to  compensate  airspeed  during  extreme 
shears,  is  proposed.  Such  a configuration  has  been  shown  to  yield  signi- 
ficant improvements  during  knife-edge  shears. 
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The  altitude  command  ( h)  function  of  range  ( R)  was  established 
to  fix  in  space  a path  similar  to  that  of  the  nominal  exponential  flare,  in  the 
presence  of  zero  winds.  The  function  was  derived  as  follows; 


Let; 


Constraints;  at  h = 60'  dh 

R = 2250'  ’ dr  ^ 

at  h = O'  . ^ _ 

R = 200' ’ dr 


2 ft  * 

g'  = . 0467  ( rad) 

0.5  (deg.)  = .0087  (rad) 


( 1)  0 = a(  200)  2 + b(  200)  + c 

( 2)  2a(  2250'  ) + b = . 0467 

( 3)  2a(  200'  ) + b = . 0087 


subtracting  ( 3)  from  ( 2) ; 


2a  ( 2050)  = .0467  - .0087  = .0380 

.*.  a = 9.3x1 0"^ 
b = .0062 


from  equation  ( 1 ) ; 

c = -1.6 


use 

A 

h = 9.3 

/ R \ ^ + 6.  2 , 

( R > 

- 1.6 

looo;  1 

^ 1000  ^ 

This  function  has  been  plotted  in  Figure  17. 
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